Nondiffusive heat transfer in muscle tissue. Preliminary results 
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We present preliminary experimental data that enable us to suggest that heat transfer in cellular 
tissue under local strong heating is a more complex phenomenon than a simple heat diffusion. 
Namely, we demonstrate that under local strong heating of a muscle tissue heat transfer in it 
exhibits substantial anisotropy unexplained in the context of the standard diffusion model. The 
observed temperature dynamics is also characterized by nonlinear behavior as well as by a certain 
repeat reversibility. The latter means that the time variations in the temperature of a cellular tissue 
undergoing repeated acts of heating go in the same way at least approximately. We explain the 
observed anomalous properties of heat transfer by suggesting the flow of the interstitial liquid to 
appear due to nonuniform heating which, in turn, affects the heat transfer. A possible mechanism 
responsible for this effect is discussed. 



I. INTRODUCTION 

During the last two decades the description of heat 
transfer in living tissue (biohcat transfer) and mathe- 
matical modeling of temperature distribution in it under 
local strong heating have attracted attention due to sev- 
eral reasons. First, this problem is an important part 
of planning and optimizing hyperthermia treatment and 
thcrmotherapy of tumors where the precise control over 
the tissue temperature is essential (see, e.g., Jl|, g, [§). 
Second, what is interesting from the standpoint of fun- 
damental physics, living tissue is a complex medium in 
which highly branching blood vessels form a fractal net- 
work of fast heat transport j| , endowing heat propaga- 
tion in it with anomalous properties g]. Third, living 
tissue is an active medium, in particular, it responds to 
strong heating, which causes the blood perfusion rate to 
grow locally by tenfold (6). In this case one meets also 
the basic problem of how a natural hierarchical system 
comprising a huge amount of elements can respond prop- 
erly to changes in its state when none of the elements 
possesses the total information required to govern the 
system as a whole [0, |). 

By now a number of significant steps towards devel- 
oping the bioheat transfer theory have been made. In 
particular, Chen & Holmes ||] actually stated the mi- 
croscopic description of the problem. Baish et al. pfj| , 
Lagendijk et al. [[J, Roomer et al. ||, 0, and 
Gafiychuk & Lubashevsky || studied the effect of the 
vessel discreteness on the temperature field. Wcinbaum 
et al. |l5|, [l6| and Lubashevsky & Gafiychuk ||] took 
into account sophisticated details of heat interaction be- 
tween different vessels, namely, the counter-current ef- 
fect, Chato |L7| and Weinbaum & Jiji |L5]], that renor- 
malizes the kinetic coefficients of heat propagation in the 
tissue. Lubashevsky & Gafiychuk || [| proposed a hi- 
erarchical model for the tissue response to local heating. 
There is also an approach by Lagendijk et al. (see, e.g., 
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|i8| and references therein) that dealing with individual 
vessels tackles the biohcat transfer problem numerically. 

All these models and approaches are based on the as- 
sumption that the cellular tissue is a passive medium 
and heat propagation in it obeys Fick's law, i.e. the 
heat flux Q in the cellular tissue is proportional to the 
temperature gradient VT. Namely, Q = kVT and the 
thermal conductivity k is isotropic or, at least, is slightly 
anisotropic. Indeed, according to the experimental data 
summarized in [ fl9| , for example, the value of the tem- 
perature diffusivity D = n/{cp) (c and p are the spe- 
cific heat and density of the tissue) is approximately the 
same for different tissues, D ■ 10 3 « 1.8, 1.3, 1.4, 1.5, 
1.2 cm 2 /s for human muscle, kidney, spleen, liver, and 
blood, respectively. We note that models developed in 
food engineering also regard the cellular tissue as a pas- 
sive medium (see, e.g., Q). This assumption is justi- 
fied by a number of self-consistent results obtained in 
a variety of experimental works devoted to measuring 
thermal characteristics of perfused and nonperfused bio- 
logical tissues (see, e.g., |19| and references therein, and, 
also, (2l], |22|, ^3|, ^4|, ^5|). However, measuring thermal 
characteristics of living tissue is a difficult task because of 
its complex structure. So, interpretation of the obtained 
experimental data and verification of the corresponding 
mathematical models are related problems which are far 
from being completed (see, e.g., |p6[ p7f), concerning es- 
pecially strong heating. Nevertheless, up to now all the 
peculiarities of heat transfer in living tissue are assumed 
to be due to the convective transport with blood. 

During hyperthermia treatment the tissue temperature 
grows up to the upper boundary of the tissue surviving, 
i.e. up to 44 — 46 °C, moreover, during a thermother- 
apy course based on the thermal tissue coagulation the 
temperature can attain values about 100 °C. So the ques- 
tion of whether the cellular tissue undergoing so strong 
heating can be described in the framework of the passive 
medium model is not evident. Therefore the thermal 
properties of cellular tissue under local strong heating 
deserve individual investigations to clarify if the physi- 
cal background of the bioheat transfer theory has to be 
modified. In particular, the fact that heat propagation 
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FIG. 1: Scheme of measuring the tissue temperature inside 
tissue bulk under local heating caused by current running 
through the resistance. 



in cellular tissue can exhibit anomalous properties was 
noted in |§. 

In this paper we present preliminary experimental data 
that have enabled us to doubt whether cellular tissue, at 
least some of its types, is in fact a passive medium and 
heat propagation in it does obey Fick's law. 



II. TEMPERATURE DYNAMICS IN THE 
TISSUE BULK UNDER STRONG LOCAL 
HEATING 

First, we studied time variations of the temperature 
inside a fragment of cow haunch muscle not undergone 
freeze after death and cut out so that muscle fibers be 
parallel to its surface (Fig. |l|). A small resistance r of 
size 1x2 mm was embedded into the tissue at a depth 
of 3 cm through this free surface oriented horizontally. 
Passing a current / through it we got a local heat source 
of controlled power. Thermocouples measuring the tem- 
perature were located at a distance p of 1 cm from the 
resistance along the fibers and in the transverse direction. 

We studied the heating dynamics for various values 
of the initial tissue temperature from 0°C to 30 °C and 
various specimens of muscle, the results were practically 
identical. Fig. demonstrates typical time variations of 
the tissue temperature obtained for two values of the heat 
source power equal to 1.2 W (/ = 300 mA, r = 13 0) 
and 2.7 W (I = 300 mA, r = 30 O). The time origin, 
t = 0, is the beginning of heating and the vertical dashed 
line, t ~ 20min, represents the heating tune-out. All the 
shown curves were obtained for tissue regions undergone 
heating for the first time because under strong heating 
irreversible transformations of the living tissue structure 
can occur. It should be noted that the curves correspond- 
ing to heating power of 1.2 W were obtained for the tissue 



□ ) source power m \ source power 

O / of 2.7 W -•- J of 1.2 W 




■ i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i 

5 10 15 20 25 30 35 

Time (min) 



FIG. 2: Time variations of the tissue temperature at points 
shown in Fig. [j] for two values of heating power. The data have 
been obtained for tissue regions heated for the first time. The 
curves corresponding to heating power of 1.2 W are shifted 
down by 1.4 °C to make the comparison more easy. 



with the initial temperature Ti n j t = 20.2 °C, so in Fig. ^ 
they are shifted down by 1.4 °C to make the comparison 
more easy. 

As seen in Fig. ^ the time variations of the tissue tem- 
perature along the fibers and in the transverse direction 
exhibit substantially different features. In fact, for the 
longitudinal direction the characteristic time 77 of the 
temperature variations can be estimated as 77 ~ 6 — 8 min 
from the e-fold decrease after the heating tune-out. The 
same estimate will be get if we apply to the point 
where the curve describing the temperature increase un- 
der 1.2 W heating transfers from the concave to convex 
form. For the temperature increase under 2.7 W heating 
the latter time seems to be even shorter. The analo- 
gous curves describing the temperature variations in the 
transverse direction are characterized by a substantially 
longer time scale 77. In particular, the temperature in- 
crease curves remain concave up to the heating tune-out. 
So, if we related the time scales 77,4 to the temperature 
diffusivities Di t along the fibers and transversely them 
with the expression 77.4 ~ p 2 /Di it then we would have 
been to set that Di/D t > 3. From our point of view 
the latter estimate contradicts the physical properties of 
heat transfer in cellular tissue because it does not contain 
any regions where the heat conductivity takes so different 
values that could give rise to this estimate. 

Under strong heating the cellular tissue coagulates and 
changes its properties, however, on the time scales un- 
der consideration it is essential only for temperatures 
exceeding 50 — 60 °C p8]| . Therefore, this effect seems 
not to be substantial enough to play a key role in the 
heat propagation through the cellular tissue in the given 
case. Indeed, if we again assume the tissue temperature 
to be governed by the diffusion equation with, may be, 
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FIG. 3: Time variations of the tissue temperature at points 
shown in Fig. |l]for the cellular tissue undergone two- fold heat- 
ing of power of 1.2 W. 



FIG. 4: Scheme of testing heat propagation through a tissue 
layer bounded by wall permeable (upper figures) and imper- 
meable (lower figures) to the interstitial liquid. 



an anisotropic diffusivity, then the stationary tempera- 
ture field T caused by a <5-source will be scaled up as 
(T — T init ) — > (T — T init )A when all the spatial scales 
are reduced by the factor A, i.e. p — > p/X. So, under 
2.7 W heating the tissue temperature should attain values 
about 60 °C at points being at a distance of 0.5 cm from 
the heat source along the fibers and at a substantially 
shorter distance in the transverse direction. Thus the 
region where the cellular tissue has to coagulate during 
the heating is of the relative volume not exceeding 10% 
with respect to the lxlxl cm 3 region actually tested in 
the experiment. Under 1.2 W heating the region of the 
coagulated tissue is still less and even seems not to form 
at all because of finite dimensions of the resistance. 

Fig. |2| also shows that heat propagation in cellular tis- 
sue under strong heating must be nonlinear. At least, 
the change of the heating power from 1.2 W to 2.7W, i.e. 
by 2.3 times, corresponds to the temperature increase 
by 2.7 times at the points distant from the heat source 
along the fibers and only by 1.4 times at the points in 
the transverse direction. 

Fig. |^ demonstrates our test of the repeat irreversibil- 
ity of the heating process showing if irreversible transfor- 
mations of the tissue structure caused by strong heating 
are essential in the case under consideration. Namely, 
to answer the question of whether heat transfer in cellu- 
lar tissue exhibiting such anomalous properties can go in 
the same way under repeated acts of heating we tested 
the temperature dynamics for the source power of 1.2 W. 
The obtained curves present time variations of the tissue 
temperature at the points shown in Fig. ^ during two 
successive processes of heating of the same tissue region 
that are separated in time by 40 min. We see not too sig- 
nificant change in the temperature variations only at the 
point distant from the heat source along the fibers that 
occurs at the stage of direct heating. During the repeated 
heating the temperature increase is characterized solely 



by a faster growth rate. So the mechanism responsible 
for these anomalous properties of heat transfer seems to 
be of reversible nature, at least, at first approximation. 

We assume that these heat transfer anomalies are due 
to motion of the interstitial liquid induced by heating, 
which is demonstrated in the next section. 



III. TEMPERATURE DYNAMICS IN A TISSUE 
LAYER UNDER STRONG SURFACE HEATING 

Following actually a standard way of measuring the tis- 
sue heat conduction we studied heat propagation through 
a tissue layer of thickness of 2 cm (Fig. [|). One of its 
surfaces was kept up at high temperatures about 100 °C 
whereas the temperature growth at the opposite surface 
was recorded by a thermocouple. The tissue layer was 
oriented vertically to avoid the gravitation effect. We 
studied four arrangements distinguished by the fiber ori- 
entation with respect to the walls bounding the layer 
and by the wall permeability to the interstitial liquid. In 
one case we used the impermeable walls to suppress the 
induced motion of the interstitial liquid, in the other we 
applied a grid to make it as free as possible. The distance 
between the walls was controlled by a certain small fixed 
force pressing the system so that there were no air-gaps 
between the tissue and the walls. 

It should be noted that in the given situation heat 
propagation through the cellular tissue differs substan- 
tially from the case considered in the previous section. In 
fact, when the thickness of the tissue layer is much less 
than the other dimensions and its surface is heated uni- 
formly, the temperature field can practically vary solely 
along the direction perpendicular to the layer surface. 
Under such conditions the interstitial liquid flow can 
appear only in the same direction and inevitably has 
to be affected by the properties of the opposite layer 



4 





90- 




og 


80 -j 




70 




<u 
H 
P 


60 i 




tperat 


50- 




Teir 


40 -i 
30 -i 






20 -j 





I 1 1 1 I 1 1 1 I 1 1 1 I 1 



' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' 1 ' I ' ' ' I ' ' ' 1 1 

0123456789 10 



Time (min) 



heat source 

1 cm drop of aniline dye 




tissue 



FIG. 6: Scheme of the experiment demonstrating the mass 
propagation induced by nonuniform heating of the cellular 
tissue. 



FIG. 5: Time variations of the temperature at the boundary 
of the tissue layer opposite to the heated one (Fig. Q). To make 
the relation between the specific data and the corresponding 
tissue-wall configuration clear the latter is labeled in Fig. W 
with the same symbol. 



surface. For local heating of the cellular tissue at in- 
ternal points distant far from its surface the structure 
anisotropy should give rise to the interstitial liquid flow 
of partly swirling form. In this case the effect of the 
surface properties seems to be weakened. 

Fig. H exhibits the obtained results. When the tissue 
layer is bounded by the impermeable walls and the inter- 
stitial liquid has no way to move the evolution of tissue 
temperature is to be governed by the conventional heat 
diffusion and to depend weakly on the fiber orientation. 
Exactly such a behavior was observed. Then assuming 
the tissue temperature to obey the diffusion equation we 
can directly estimate the corresponding temperature dif- 
fusivity D. Namely, in the given case the tissue tempera- 
ture at one boundary is fixed and equal approximately to 
100 °C and at the other we may impose the zero flux con- 
dition on the temperature field. The initial temperature 
is 20 °C and according to Fig. [| the tissue temperature 
attains values about 50 °C at latter boundary in 10 min. 
Then fitting the solution of the diffusion equation sub- 
ject to the appropriate conditions to these data we get 
the estimate D ~ 2 x 10 -3 cm 2 /s. The latter corresponds 
well to the available data for the tissue temperature dif- 
fusivity within the accuracy of our measurements. 

In the tissue layer bounded by the grid heat transfer 
exhibits another behaviour. When the muscle fibers were 
parallel to the layer boundaries the temperature growth 
turned out to be of the same form as for the impermeable 
walls. By contrast, when the tissue fibers were perpendic- 
ular to the boundaries the temperature growth deviated 
substantially from the previous dependence (Fig. |5|) . We 
can explain this effect only admitting an essential flux of 
the interstitial liquid to appear due to the nonuniform 
heating of the cellular tissue. It is likely that the fiber 
arrangement gives rise to long channels in the muscle 



tissue running along the fibers. So, when there are no 
substantial obstacles to motion of the interstitial liquid 
nonuniform heating of the tissue is able to induce the 
interstitial liquid flux affecting remarkably heat trans- 
fer. To partly justify the given assumption we note that 
for the latter fiber orientation the blood efflux through 
the grid was visually strong in comparison with the other 
cases when the measured tissue temperature attained the 
values about 33 °C. 

In order to make this phenomenon more evident the 
next section presents its certain visualization. 

IV. MASS PROPAGATION CAUSED BY 
NONUNIFORM HEATING 

If nonuniform heating of cellular tissue causes the in- 
terstitial liquid to flow in it affecting substantially heat 
propagation itself then the interstitial liquid flux can be 
directly visualized through the induced propagation of a 
certain substance dying the tissue. Indeed, diffusion co- 
efficients of simple molecules in cellular tissue (of order 
10~ 5 cm 2 /s or even less) are much less than the temper- 
ature diffusivity, so mass propagation in cellular tissue 
is to be essentially accelerated by this interstitial liquid 
flux. 

Figure ^| illustrates the experiment showing the effect 
of nonuniform heating on propagation of aniline dye in 
the cellular tissue. We put an aniline dye drop in a lunula 
on the horizontal tissue surface and pressed down the tis- 
sue with a glass plate to make the tissue surface plane 
and, thereby, to exclude the influence of the surface re- 
lief. A metallic disc of diameter of 1 cm was embedded 
into the glass plate placed so that this disc be right over 
the drop. An ohmic resistance was attached to the ex- 
ternal disc surface and passing a current through it we 
locally heated the tissue. The colour patterns formed by 
the aniline propagation are represented in Fig. 0. Fig- 
ure (a) exhibits the colour pattern appeared in the tissue 
without heating after a lapse of 8 hours. As is seen, 
it is roughly of symmetrical form and the muscle fibers 
have no substantial effect on the aniline propagation in 
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FIG. 7: Colour patterns formed by aniline dye spreading over 
the tissue from the initial small drop (Fig. []) without heating 
(a) and under strong local heating (b). (In the this figure the 
initial drop of aniline dye located under the heating source is 
invisible.) 

this case. The opposite situation is shown in figure (&) 
for the aniline propagation under a strong local heating. 
Already after 15 min the colour pattern attained the pre- 
vious one in size and got a substantially anisotropic form, 
the pattern dimension along the fibers is more that twice 
its mean size in the transverse direction. Moreover, in 
this figure it is clearly seen as aniline dry propagates in- 
side small grooves on the tissue surface whereas without 
heating such grooves practically have no effect. 



V. DISCUSSION AND HYPOTHESES 

The presented preliminary experimental data enable 
us to suggest that heat transfer in cellular tissue un- 
der local strong heating is a more complex phenomenon 
than a simple heat diffusion. On one hand, the ob- 
served strong anisotropy in the heat propagation along 
the muscle fibers and in the transverse direction cannot 



be explained keeping in mind only the heat conduction 
through the cellular tissue. Indeed, otherwise, we should 
assume that the temperature diffusivity along the fibers, 
Di, is fully three times greater than the corresponding 
value in the transverse direction, Dt (Di > Dt), whereas 
the cellular tissue contains no elements differing from one 
another so remarkably in the heat characteristics. On 
the other hand, the fiber arrangement endows the mus- 
cle tissue with a significant anisotropy for motion of the 
interstitial liquid through the tissue as through a porous 
medium. Therefore, if the nonuniform temperature field 
can cause the interstitial liquid to flow then heat propaga- 
tion in the muscle tissue will be able to exhibit substantial 
anisotropy in properties. We think that in the presented 
experiments exactly this phenomenon, i.e., the intersti- 
tial liquid flux induced by nonuniform heating and, in 
turn, its effect on heat transfer was observed. Obviously, 
it should be nonlinear, which corresponds to the data 
presented in Fig. |[ 

The experimental data shown in Fig. |^ enable us to as- 
sume that heat transfer in cellular tissue under such con- 
ditions is characterized by the repeat reversibility. The 
latter means that the time variations in the temperature 
of a cellular tissue undergoing repeated acts of heating 
go approximately in the same way. So a possible micro- 
scopic mechanism responsible for the induced flux of the 
interstitial liquid is also of the reversible nature, at least, 
it is not due to the tissue thermal coagulation or another 
thermal damage of the tissue. For example, it could be 
implemented through the liquid redistribution between 
the tissue cells and the intercellular space due to the 
change in the osmotic equilibrium induced by heating. 
In the given case an additional amount of the interstitial 
liquid appearing in the heated region has to give rise to 
its flow from this region. 

We point out again that the interstitial liquid flux in- 
duced by heating seems to be the integral feature of heat 
transfer in cellular tissue. Therefore, only under special 
conditions heat propagation in cellular tissue is reduced 
solely to diffusion. Moreover, mass and heat propaga- 
tion have to be related with each other through this flux, 
which is illustrated in Figs. || and 0. 

Concluding the paper we would like to note that we 
observed anomalous properties of heat transfer in muscle 
tissue where the fibers endow the tissue with well devel- 
oped anisotropy. There are organs, for example, kidney 
and liver whose cellular tissue does not possess signifi- 
cant anisotropy on the average regarding such an organ 
as a whole. For rigorously isotropic media the effect of 
a similar liquid flux induced by heating is likely to be 
depressed. However, on smaller scales these organs are 
certainly to exhibit local anisotropy with random orien- 
tation in space. Therefore, local strong heating again has 
to induce the interstitial liquid flux that, may be, changes 
its direction randomly in different regions of a given or- 
gan. In this case heat transfer is also to exhibit anoma- 
lous properties, which cannot be described, in general, 
by the conventional diffusion equation (see, e.g., pf][). 
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